Riemannian thermo-statistics geometry 



L. Velazquez 

Abstract. It is developed a Riemannian reformulation of classical statistical me- 
chanics for systems in thermodynamic equilibrium, which arises as a natural ex- 
tension of Ruppeiner geometry of thermodynamics. The present proposal leads 
to interpret entropy S g (I\8) and all its associated thermo-statistical quantities as 
purely geometric notions derived from the Riemannian structure on the manifold of 
macroscopic observables M.g (existence of a distance ds 2 = gij(I\6)dPdP between 
macroscopic configurations I and I + dl). Moreover, the concept of statistical cur- 
vature scalar R{I\9) arises as an invariant measure to characterize the existence 
of an irreducible statistical dependence among the macroscopic observables / for a 
given value of control parameters 6. This feature evidences a certain analogy with 
Einstein General Relativity, where the spacetime curvature i?(r, t) distinguishes 
the geometric nature of gravitation and the reducible character inertial forces with 
an appropriate selection of the reference frame. 
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1. Introduction 



Let be a classical equilibrium distribution function: 
dp(I\9) = p(l\0)dl, 



(1) 



where the stochastic variables I = {P} denote the relevant macroscopic observables 
of the system under analysis, while 9 = {9 a } the underlying control parameters of 
a given equilibrium situation. It is possible to introduce a distance notion between 
two close equilibrium situations 9 and 9 + d9: 



The existence of this type of Riemannian formulation was pioneering proposed by 
Rao j2], which is referred to as inference geometry in the literature |3j. Alternately, 
it is also possible to introduce a distance notion between two close macroscopic 
configurations I and I + dl for a given value of control parameter 9: 



starting from the same distribution function ([I]). This latter geometric 
characterization appears as a suitable extension of Ruppeiner geometry of 
thermodynamics [I] to the framework of equilibrium classical statistical mechanics. 
For this reason, the present proposal shall be hereafter referred to as Riemannian 
thermo- statistics geometry. The main interest of this work is to present a systematic 
development of the most relevant features of this geometric formulation starting from 
an axiomatic perspective. Such a procedure allows to show that fundamental concepts 
of statistical mechanics can be suitably rephrased in terms of geometry notions, which 
provides a general framework to apply powerful methods of Riemannian geometry 
for the analysis of properties of thermodynamical systems. 



ds 2 F = g a0 (9)d9 a d9 
in terms of the so-called Fisher's information matrix p]: 



(2) 




(3) 



ds 2 = gij (I\9)dPdP 



(4) 
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2. Riemannian thermo-statistics geometry 

Let us denote by Ai the abstract manifold composed of all admissible values of 
macroscopic observables /, as well as by V the abstract manifold composed of all 
admissible values of control parameters 9. Besides, let us also denote by Atg the 
sub-manifold of Ai composed of all values of macroscopic observables / that are 
accessible for a given value 9. In general, it is possible to consider two different types 
of coordinate reparametrizations: (1) the coordinate reparametrizations Q(I) : Hi — > 
TZq of the manifold of macroscopic observables Ate, as well as (2) the coordinate 
reparametrizations u{9) : He — > H v of the manifold of control parameters V . 

2.1. Postulates of thermo-statistics geometry 

Axiom 1 The manifold of the system macroscopic observables Aie possesses a 
Riemannian structure, that is, it is provided of a metric tensor (I\9) and a 
torsionless covariant differentiation D; L that obeys the following constraints: 



Definition 1 The Riemannian structure on the manifold of macroscopic observables 
Aie allows to introduce the invariant volume element as follows: 



where \gij (I\9)\ denotes the absolute value of the metric tensor determinant. 

Axiom 2 There exist a differentiate scalar function S g (I\9) defined on the manifold 
M. e , hereafter referred to as the scalar entropy, whose knowledge determines the 
equilibrium distribution function dp(I\9) of the macroscopic observables I e Aig as 
follows: 



Definition 2 Let us consider an arbitrary curve given in parametric form I(t) G Ai$ 
with fixed extreme points I(ti) = P and 1(^2) = Q- Adopting the following notation: 



D k g l3 (I\6) = 0. 



(5) 




(6) 



dp (I\9) = exp [S g (I\9)} du(I\6). 



(7) 



P = 



dP(t) 



(8) 



dt 
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the length As of this curve can be expressed as: 

<-t 2 



As= / yJ gij [I(t)\e}P(t)P(t)dt (9) 

Jt\ 

Definition 3 The curve I(t) G A4q exhibits an unitary affine parametrization 

when its parameter t satisfies the following constraint: 

g tJ (i\e)i l (t)P(t) = i. (io) 

Definition 4 A geodesic is the curve I g (t) with minimal length (Tl\) between two 
fixed arbitrary points (P, Q) G M.q. Moreover, the distance D(P, Q\8) between these 
two points (P,Q) is given by the length of its associated geodesic I g (t): 

D(P, Q\9) = I" y/ 9ij [I g (t)\9]i g (t)P g (t)dt. (11) 

Definition 5 Considering a differentiable curve I(t) G M.q with an unitary affine 
parametrization, the entropy production along this curve $(t) is given by: 

(JjL 

Axiom 3 The length As of any interval (t\, t 2 ) of an arbitrary geodesic I g (t) G M.q 
with an unitary affine parametrization is given by the variation of its entropy 
production A$(t) with opposite sight: 

As = -A$(t) = $(ti) - $(t 2 ). (13) 

Axiom 4 The probability density p(I\0) associated with distribution function |7p 
vanishes with its first partial derivatives for any point on the boundary dM. e of the 
manifold Me- 

2.2. Interpretations and fundamental consequences 

Commentary. Axiom 1 allows to precise the Riemannian structure of the 
manifold Aie starting from the knowledge of the metric tensor gij(I\8), specifically, 
the covariant differentiation Di and the curvature tensor Rijki(I\6). As discussed 
elsewhere [7j, Eq.((5]) is an strong constraint of Riemannian geometry that determines 
the affine connections employed to introduce the covariant differentiation D i: 
specifically, the so-called Levi-Civita connection: 

lT \a\ _ „km l ( d 9im . dg jm dgij 
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The knowledge of the affine connections allows the introduction of the curvature 
tensor R\- k = R\- k (I\9) of the manifold Me'- 

8 8 

T)l W _ p/ _i_ p£ pm p/ prra /-i r\ 

ijk 8T ^ 8l3 ^ im jk jm iki V / 

which is also determined from the knowledge of the metric tensor gij(I\8) and its first 
and second partial derivatives. Using the curvature tensor R\- k (I\9), it is possible to 
obtain its fourth-rank covariant form Rij k i(I\9) = gi m (I\9)R™ k (I\9): 

f> - 1 ( d2gu d2 ^ k d 2 g j7 d 2 g ik \ 
ijkl 2 \8V8I k 8P8I 1 dPdP 8V8I 1 ) 

+ 9mn (r™r™ fc — r^r™ fc ) , (16) 

the Ricci curvature tensor Rij(I\6): 

Rtiim = R k kl3 {I\6) (17) 
as well as the curvature scalar R(I\6): 

R{I\0) = g^mRUm = g ij (I\e)9 kl mRkiji(I\0)- (18) 

The curvature scalar R(I\6) has a paramount relevance in Riemannian geometry [7j 
because of it is the only invariant derived from the first and second partial derivatives 
of the metric tensor gij(I\6). ■ 

Commentary. Axiom 2 is a covariant redefinition of Einstein postulate of classical 
fluctuation theory [5]: 

dp EP (I\9) = Aexp [S(I\9)] dl, (19) 

which rephrases Boltzmann entropy S = log W to assign relative probabilities from 
the entropy S(I\9). However, expression ([19]) has the disadvantage that the entropy 
S(I\9), commonly referred to as the coarsed grained entropy, does not correspond to 
a scalar function within a geometric theory. In fact, coarsed grained entropy S(I\9) 
behaves under coordinate reparametrizations : Hi —> He as: 

8Q 



S(Q\9) = S(I\9)- \og 



(20) 



81 

The non-scalar character of the entropy S (I\9) is inconsistent with the physical idea 
that thermodynamic entropy is a state function, which should not depend on the 
particular coordinate representation employed to describe the macroscopic properties 
of a thermodynamic system. One can avoid this mathematical inconsistence replacing 
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the usual volume element dl by the invariant volume element (j6]). Thus, the scalar 
character of the entropy S g (I\9): 

S g (Q\9) = S g (I\9) (21) 

and the covariance of the metric tensor (I\9): 

dl m dl n 

= ^ W 9mn{I\e) (22) 

guarantee the invariance of the equilibrium distribution function (j7J) under coordinate 
reparametrizations 6(7) : IZj — >■ TZq. ■ 

Remark 1 (Alternative Axiom 2) Starting from the knowledge of probability 
density p(I\9) of the distribution function (T7]j and the existence of an everywhere 
non-vanishing metric determinant \gij(I\9)\, it is possible to introduce the scalar 
entropy S g (I\9) as follows: 

S g (I\9) = log . =. (23) 

Commentary. Axiom 3 states a direct relation between the distance notion (j3J) 
and the entropy production (1531 . which relates the metric tensor gij(I\9) and the 
scalar entropy S g (I\9). This postulate can be regarded as a generalization of the 
Ruppeiner postulate [I]: 

a (1\9\- 92Sa m (24) 

which provides a metric tensor for thermodynamics, with I being the point with 
global maximum entropy. ■ 

Theorem 1 The scalar entropy S g (I\9) is locally concave everywhere and the metric 
tensor gij(I\9) is positive definite on the manifold Me- Moreover, the metric tensor 
gij(I\9) can be identified with the covariant Hessian f-tij(I\9) of the scalar entropy 
S g (I\9) with opposite sign: 

g l3 (I\9) = -7Uj{I\e) = -DiDjSg (I\9) . (25) 

Proof. The searching of the curve with minimal length (Till between two arbitrary 
points (P, Q) is a variational problem that leads to the following ordinary differential 
equations [7]: 

i k g (t)D k r g (t) = r g (t) + r mn [i g (t)\9] i™{t)i n g {t) = o, (26) 
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which describes the geodesic I g (t) with an unitary affine parametrization. Eqs. (!83|) 
and (fT3j) can be rephrased as follows: 

at/ \ d$(s) d 2 S„ , . „, 

As = _A*(.) - -Ai = = -!. (27) 

Taking into account the geodesic differential equations (|26|) . constraint ( 127]) can be 
rewritten as: 

d 2 s g = j k ds g fifj d 2 s g = , ifj r d 2 s 3 fe as g ] 

ds 2 5 99 dPdP 9 9 \ dPdP lJ dl k J ' 1 ; 

where it is possible to identify the covariant entropy Hessian Tiif 

«« = W 3 = |4 - r* §. (29) 

Eqs. (l28|) - (|29|) can be combined with constraint (|T0|) to obtain the following 
expression: 

fa* + = o- (so) 

Its covariant character leads to Eq. fl25l) . The concave behavior of the scalar entropy 
S g (I\9) and the positive definition of the metric tensor gij(I\8) are two direct 
consequences of Eq. (f2"Tj) . ■ 

Corollary 1 The metric tensor g^ (I\9) can be obtained from a given scalar entropy 
S g (I\9) through the following set of covariant partial differential equations: 

= d 2 S g 1 km ( dgirn Ogjrn _ %A OS^ 

yiJ " dPdP 2 g \dP dP dl m J dP' 1 ' 

The admissible solutions derived from the nonlinear problem (3l\) should be 
everywhere finite and differentiate, including also on boundary of the manifold Aig. 

Corollary 2 Ruppeiner geometry of thermodynamics /^/ is obtained from as 
a particular case of thermo- statistics geometry restricting to the gaussian 
approximation of the distribution function (0): 



dp (I\9) ~ J\f exp 



-\ gij {l\9) 5P5P 



\ 



9* {AO) 



2tt 



dl, (32) 



where 51 = I— I andM = exp \S g (-f |$)] ■ Moreover, the normalization of distribution 
function [3^) leads to the following estimation for the maximum scalar entropy 
Sg{l\9): 

S g (I\9) ~ 0. (33) 
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Proof. Eq. (l25|) drops to Ruppeiner definition of thermodynamic metric tensor (J2] 
restricting to the point / with of the global maximum scalar entropy S g (I\6). For 
a sufficiently large thermodynamic short-range interacting system, it is possible to 
assume the following scaling dependencies with the characteristic size parameter 
A, S g ~ A, / ~ A =>- gij ~ 1/A and the correlation functions (SPSP) ~ A, 
where the role of A can be performed by the system volume V, the total mass 
M, number of constituents N, etc. Gaussian distribution fl32l) represents the 
lowest approximation level that accounts for the thermodynamic fluctuations of the 
macroscopic observables / in the power expansion of the scalar entropy S g (I\6) and 
the logarithm of the factor yj\gij(l\6)\. ■ 

Commentary. Axiom 4 talks about the asymptotic behavior of the equilibrium 
distribution function ([7]) for any point h on the boundary OM.g\ 



lim p{I\ff) = Jim ^p(I\e) 
I-*h i->h or 



0. 



(34) 



Such conditions play a fundamental role in the character of stationary points 
(maxima and minima) of the scalar entropy S g (I\6) as well as the fluctuation 
theorems derived on the Riemannian structure of the manifold Aig. Due to their own 
importance, such fluctuation theorems shall be discussed in a forthcoming paper. ■ 

Remark 2 The boundary conditions (34\ ) are independent from the admissible 
coordinate representation TZj of the manifold Aig. 

Proof. This remark follows as a direct consequence of the transformation rule of 
the probability density: 

-i 



P (e\9) = P (i\9) 



01 



as well as the ones associated with its derivatives: 

d P {e\o) _ dp f d P (i\o) d 



dQ 



01 



oe 



01 



(35) 



(36) 



under a coordinate reparametrization 0(1) : TZi —> TZq whose Jacobian \8Q/dI\ be 
everywhere finite and differentiable. ■ 
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3. Gaussian and spherical representations 

3.1. Gaussian Planck potential V g {6) 

Definition 6 The covariant form of the gradiental generalized forces ipi (I\9) 
are defined from the scalar entropy S g (I\9) as follows: 

V< (I\6) = -DiSg (I\9) = -dS g (I\9) /dP. (37) 

Using the metric tensor g l i (1 \9) , it is possible to obtain its contravariant counterpart 
^(I\9): 

^(I\9)= g v(I\9)^(I\9), (38) 
as well as its the square norm ip 2 = ip 2 (I\9): 

^ 2 (I\9) = ^(I\9)^(I\9). (39) 

Theorem 2 The scalar entropy S g (I\9) can be expressed in terms of the square norm 
of the gradiental generalized forces as follows: 

S g {I\9) = V 9 {9)- l -^{I\9\ (40) 

where V g {9) is a certain function on control parameters 9 hereafter referred to as the 
gaussian Planck potential. 



Proof. Let us introduce the scalar function V g (I\9): 

V g (I\9) = S g (I\6) + l -gV{I\9)Ul\0)Um- (41) 
It is easy to verify that its covariant derivatives: 

D k V g {I\9) = D k S g {I\9) + i {^{I\9)^ j {I\9)D k g i '{I\9)+ (42) 
+g ij {I\9) Ml\9)D k ^{I\9) +ij Jj {I\9)D k ^ i {I\9)]} 

vanish as direct consequences of the metric tensor properties (jSJ) and (J2"5j) . as well as 
the definition f|3T|) of the gradiental generalized forces ipi(I\9). Since the covariant 
derivatives of any scalar function are given by the usual partial derivatives: 

D kP a m = £k'P*(W = 0, (43) 
the scalar function V g (I\9) only depends on the control parameters: 

p g (I\0) = VM (44) 
This last result leads to Eq. lHUI) . ■ 
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Corollary 3 The value of scalar entropy S g (I\9) at all its extreme points derived 
from the stationary condition: 

^ 2 {I\9) = (45) 

is exactly given by the gaussian Planck potential V g (9). 

Corollary 4 The equilibrium distribution function |7]) admits the following 
gaussian representation: 

1 



dfi(I\9). (46) 

Here, the factor Z g {9) is related to the gaussian Planck potential as follows: 

V 9 (9) = - log Z g (9), (47) 
which shall be hereafter referred to as the gaussian partition function. 

3.2. Maximum and completeness theorems 

Theorem 3 The scalar entropy S g (I\9) exhibits a unique stationary point I, which 
corresponds to its global maximum. 

Proof. Since the metric tensor gij(I\9) defines a positive definite metric on the 
manifold A4g, the vanishing of the probability density: 



p{I\9) = exp [S g (I\9)) 



9a m 



(48) 



2tt 

on the boundary dAig of the manifold A4q for every admissible coordinate 
representation TZi evidences the vanishing of the scalar function: 

*(/|0) = exp[<S g (/|0)] (49) 

on the boundary dAig. Since the function ty(I\9) is nonnegative, finite and 
differentiable on the manifold A4$, the entropy S g (I9) should exhibit at least a 
stationary point where takes place the stationary condition f )45p . Since the scalar 
entropy S g (I9) is a concave function, its stationary points can only correspond to 
local maxima. Let us suppose the existence of a least two stationary points I\ and I2, 
which can always be connected with a certain geodesic I g (t). According to constraint 
(|27|) . the entropy production $(t) is a monotonous function along the curve I g (t). 
Consequently, $(t) should exhibits different values at the stationary points Ii and 
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I2, which is absurdum since the entropy production $(£) identically vanishes for any 
stationary point of the scalar entropy S g (I\9): 

9{t) = -P(t)MHtM- ( 5 °) 

Consequently, there exist only one stationary point that corresponds with the global 
maximum of the scalar entropy S g (I\8). ■ 

Theorem 4 Any hyper-surface of constant scalar entropy S g (I\9) is just the 
boundary of a n- dimensional sphere S n (I,r) C Aig centered at the point I with 
global maximum entropy, whose radio r = r(I\9) = D(I, I\9) at the point I allows to 
rephrase the scalar entropy as follows: 

S g {I\9) = -P g {9)- l -r\l\9\ (51) 

with n = dim Aig being the dimension of the manifold Aig. 

Proof. By definition, the vector field v l (I\9): 

is the unitary normal vector of the hyper-surface with constant scalar entropy S g (I\9). 
It is easy to verify that the vector field v l (I\9) obeys the geodesic equations fl26|) : 

v\l\9)D k v\I\9) = [51 - v\I\9)v k {I\9)] = 0. (53) 

Hence, v l (I\9) can be regarded as the tangent vector: 
dP n (s\e) 

-jt±L=v*[I g (s\eM (54) 

of geodesic family I g (s\e) with unitary affine parametrization centered at the point 
I with maximum scalar entropy S g (I\9), I g (s = 0|e) = 7, where the parameters e 
distinguish geodesies with different directions at the origin. The entropy production 
$(s|e) along any of these geodesies is given by the norm of the gradiental generalized 
forces with opposite sign: 

*(*|e) = -^P^ VMeM = -^VM^M ( 55 ) 
Considering Eq.( fl3|) . the norm ip(I\9) can be related to the length As of the geodesic 
connecting the point / with point J with maximum scalar entropy, that is, the 
distance D(I,I\9) between the points I and /: 

i>(I\9) = D(I,I\9). (56) 
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According to the gaussian decomposition ( 140|) . the hyper-surface with constant scalar 
entropy S g (I\6) is also the hyper-surface where the norm of gradiental generalized 
forces ip(I\0) is kept constant, that is, the boundary of a n-dimensional sphere 
S n (I,r) centered at the point I with maximum entropy. ■ 

Corollary 5 (Completeness) The knowledge of the metric tensor g^{l\G) and the 
point I with maximum scalar entropy S g (I\8) determines the equilibrium distribution 
function |7]). 

Commentary. The radio r(I\8) of the n-dimensional sphere S n (I,r) of the 
theorem |4] and the invariant volume element dp(I\8) are purely geometric notions 
derived from the knowledge of the metric tensor gij(I\6) and the point / with 
maximum scalar entropy S g (I\8). Therefore, the scalar entropy S g (I\6) and all its 
associated thermo- statistical quantities represents geometric notions derived from the 
Riemannian structure of the manifold Aig. ■ 

Corollary 6 The equilibrium distribution function ([?|) can be expressed as follows: 
dp(r, q\9) = - ^ J— exp (^-\ r ^j drdZ g (q\r, 9) (57) 
where a"E(q\r,8) is the hyper-surface element: 



dZ g (q\r,9) 



9ap(r,q\6) 



2tt 



dq. (58) 



obtained from the metric tensor g a /3{ r ,(l\@) associated with the projected 
Riemannian structure on the hyper-surface dS^ n \l\r) of the n-dimensional 
sphere S^(I\r) C M g . 

Proof. Let us consider the geodesic family I g (s;e) derived as a solution of the 
problem ( 15^1) . The quantities e = {e 1 } represent the asymptotic values of the unitary 
vector field v l (I\9) at the origin: 

e i = MmiUs; e), (59) 

s— >0 

which characterize the direction of the given geodesic I g {s; e) at the point /. These 
vectors can be parameterized as e = e(q) using the intersection point q of the geodesic 
I g {s; e) with the boundary dS n (I\r) of a n-dimensional sphere S n (I\r) C M.g of the 
theorem [4j which shall be hereafter referred to as the spherical coordinates. 
One can employ the quantities p = (r, q) to introduce a spherical representation 



CONTENTS 



13 



IZp centered at the point / with maximum scalar entropy S g (I\9). The coordinate 
reparametrization p(I) : IZi — > 1Z P is given by the geodesic family / = J*[r|e(g)]: 

dP„\r\e(q)] dP„\r\e(q)] 
v\I\0) = 9[ l W \ <(I\9) = g ^ WJ . (60) 

The new n — 1 vector fields r l a {I\9) are perpendicular to the unitary vector field 
Vi(I\9), Vi(I\9)T l a (I\9) = 0, since the vectors v l (I\9) and T l a (I\9) are respectively 
normal and tangential to the boundary dS^ n '(I\r) of the n-dimensional sphere 
S^(I\r). Consequently, the components of the metric tensor in this spherical 
coordinate representation are given by: 

g rr (r, q\9) = 1, g ra (r, q\9) = g ar (r, q\9) = 0, 

g a p(r,q\9) = gij (I\9yM\9)ri(I\9), (61) 

whose components g a p{ r -, Q\9) describe projected Riemannian structure on the 
boundary dS^ n \l\r) of the n-dimensional sphere S^ n \l\r). Eq.( l57l) is just the 
expression of the distribution function ([7j) in the spherical representation TZ p . m 

4. Statistical curvature 

Definition 7 A set of stochastic variables I exhibits a reducible statistical 

dependence when there exists at least a coordinate representation TZj of the manifold 
M.0 where the distribution function dp(I\9) can be factorized into independent 
distribution functions dp^\l^\9) = p^(I^\9)dP as follows: 

n 

dp{I\9) = Y[dp {l) {I {l) \9) (62) 
i=i 

for any I G M.q. Conversely, the set of stochastic variables I exhibits an irreducible 
statistical dependence. 

Example 1 The stochastic variables X and Y described by the distribution function: 

dp(X, Y) = A exp [-X 2 — Y 2 — XY] dXdY (63) 

are statistical dependent. However, this same distribution function can be decomposed 
into independent distributions: 



MC,v) = y ^ ex P 



-V 

2 s 



dC^=exp 

V Z7T 



1 2 



dr) (64) 
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considering the the following coordinate reparametrizations: 

X = ^=(( + V ) J Y = ^=(C-r ] ). (65) 



This fact evidences that the distribution function l[63\) exhibits a reducible statistical 
dependence. Clearly, the correlation functions: 

cov(X,Y) = ((X - (X))(Y - (Y))) (66) 

cannot provide an absolute measure of statistical dependence. 

Definition 8 A manifold A4q is flat when its geometric properties are fully 
equivalent to the n-dimensional Euclidean space W 1 . Otherwise, the manifold is said 
to be curved. 

Theorem 5 The flat character of the manifold Aie implies the existence of a 
reducible statistical dependence among macroscopic observables I , while its curved 
character implies the existence of an irreducible statistical dependence. 

Proof. According to Riemannian geometry, the flat or curved character of the 
manifold Aig is unambiguously determined by the curvature tensor Rij k i(I\9). In 
fact, the manifold Aie is flat whenever the curvature tensor R\- k {I\9) vanishes for 
any / G Aig [7]. The metric tensor gij(I\6) of an Euclidean manifold Aie can be 
written in terms of delta Kronecker matrix 8if 

9ij (I\9) = = -d 2 S 9 (I\9)/dPdP (67) 

using an appropriate coordinate representation IZi. Such a metric tensor leads to 
the following decomposition of the equilibrium distribution function dp(Q\9): 

dp(I\9) = Hexp {- [V - P(9)] 2 /2] ^= (68) 



2vr 

i 

into gaussian distributions, which evidences the reducible character of the statistical 
dependence among macroscopic observables I for the case of a flat manifold. 

To show the relationship between curved character of the manifold Ai$ and 
the irreducible statistical dependence among the macroscopic observables /, let 
us consider the spherical representation (j57]l of the distribution function (jTJ. 
The demonstration reduces to verify the irreducible coupling between the radial 
coordinate r and the spherical coordinates q for any coordinate representation lZ q of 
the boundary dS^ n \l,r). The analysis can be restricted to the asymptotic behavior 
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of the distribution functions for r small, that is, within the neighborhood of the point 
I with maximum scalar entropy S g (I\9). 

Let us consider a coordinate representation TZ X of the manifold A4e where the 
point I corresponds to the origin of the coordinate frame x = 0, and the first partial 
derivatives of the metric tensor vanish at the origin. Moreover, let us assume the 
notation A = A(pc = 0\6) to simplify the mathematical expressions. The parametric 
family of geodesic x 5 (s|e) admits an approximate solution in terms of power-series 
of the unitary affine parameter s: 



x 



Ksle) = e\q)s - -s 3 d l f i jk e\q)e k (q)e l (q) + 0(s 3 



where diT l - k is the partial derivative of the affine connection at the origin: 



Introducing the quantities Ca(l) : 



-g irn di [djgmk + d k g mj - d m g M } . 



£(?) 



de l (q) 
dq a 



(69) 



(70) 



(71) 



the components of the projected metric tensor g a p(r, q) on the boundary dS^ n \l,r) 
are given by: 

1 



g a$ (r,q) = r 2 K a p(q) - —r i R i]kl X^{q)Xf{q) + .., 



(72) 



where K a /3(q) = <7ij6a(?)£is(?) an d Rijki is the value curvature tensor ([16]) at the 
origin, while the quantities X%(q) are defined as: 

X^ q ) = e\ q )&( q )-e>( q )e a (q). (73) 

The previous expressions lead to the following approximation for the spherical 
representation of distribution function: 



dp(r, q \6) 



dg{r) 



1 - -r 2 F(q\6) + 0(r 2 ) 



dn(q). 



Here, dg(r) and dfl( q ) are two independent normalized distribution functions: 

1 



dg{r 



2f-r( t r 









dq 


2tt 


V2^ 



dr. 



(74) 

(75) 
(76) 
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while J-(q\0) is a function on the spherical coordinates q that arises as a consequence 
of non- vanishing curvature tensor Rijkf- 

JF{q\9) = R ljkl K^{q)X^{q)Xf{q). (77) 

According to Eq.(jT31). the curved character of the manifold M.q does not allow to 
decouple the radial coordinate r and the spherical coordinates q. Such a coupling is 
irreducible since the coordinate reparametrizations of the manifold A4g only affect 
spherical coordinates q in the framework of the spherical representation of the 
distribution function fl57|) . ■ 

Corollary 7 The statistical curvature tensor Rijki{I\0) allows to introduce some 
local and global invariant measures to characterize both the intrinsic curvature of the 
manifold Aig as well as the existence of an irreducible statistical dependence among 
the macroscopic variables I. They are the curvature scalar R(I\9) introduced in 
Eg. fflfy , the spherical curvature scalar H(r,q\9): 

II(r, q\9) = g^(r, q\6)R ijkl {r, q\9)X%{r, q\9)X k a l (r, q\9) (78) 

with X l J(r,q\9) being: 

X'Jir, q\9) = i/(r, g|0)^(r, q\9) - v^(r, q\9)ri(I\9), (79) 

which arises as a local measure of the coupling between the radial r and the spherical 
coordinates q in the spherical representation of the distribution function [5l\ ), and 
finally, the gaussian Planck potential V g {9) introduced in Eq.^Jlty, which arises 
as a global invariant measure of the curvature of the manifold A4$. 

Proof. As already commented, the curvature scalar is the only invariant associated 
with the first and second partial derivatives of the metric tensor gij(I\9). The 
consideration of the spherical representation of the distribution function fl57|) allows 
to introduce the normal v l (r, q\9) and tangential vectors r*(r, q\9), as well as the 
projected metric tensor g a p{r,q\6) = gtj(r, q\9)r^(r, q\9)ri{r, q\9) associated with the 
constant scalar entropy hyper-surface dS^ n \l\r). This framework leads to introduce 
the spherical curvature scalar H(r,q\9) as a direct generalization of the spherical 
function J-(q\9) of the asymptotic distribution function fITil) . 

The role of the gaussian Planck potential V g (9) as a global invariant measure of 
the curvature of the manifold A4g can be easily evidenced starting from the spherical 
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representation of the distribution function (1571) . Integrating over the spherical 
coordinates q, one obtains the following expression for the gaussian partition function: 

ZfAe) = i / + °° exp (~l r2 ) s ^ r i^ rfr ' ( 8 °) 



where H g (r\9) denotes the invariant area of the constant scalar entropy hyper-surface 
dS^(I\r). For a flat manifold, the invariant area T,fi at (r\9) is given by: 

*M<W = ^r^y (81) 

which is the area of an n-dimensional Euclidean sphere of radio r normalized by the 
factor (27r)("- 1 )/ 2 . Eq.flgUj) can be rewritten as follows: 

Z g {6) = °°ex P (~r 2 ) a(r\9)E flat (r\9)dr (82) 

introducing the geometric rate a(r\9): 

a(r\9) = X g (r\9)/X flat (r\9), (83) 

which characterizes how much differ the hyper-surface area of the sphere (I, r) 
of the manifold M. e due to its intrinsic curvature. Since the gaussian partition 
function Z g {9) = 1 for the case of the n-dimensional Euclidean space MJ 1 , a non- 
vanishing gaussian Planck potential V g (9) is a global invariant measure of the 
intrinsic curvature of the manifold Aio. ■ 



Remark 3 According to Eg. ( fffgj) and the corollary^ gaussian approximation 
employed in Ruppeiner geometry of thermodynamics [41 or % accounts for the local 
Euclidean properties of the Riemannian manifold Me at the point I with maximum 
scalar entropy S g (I\9). Consequently, the intrinsic curvature of the manifold M.e 
can be only manifested analyzing the system fluctuating behavior beyond the Gaussian 
approximation 



Theorem 6 For a sufficiently large short-range interacting thermodynamic system, 
the gaussian Planck potential V g (9) can be estimated as follows: 

V g (9) = ±R(I\9) + 0(l/A), (84) 

where R{I\9) is the statistical scalar curvature (T72j] evaluated at the point I with 
maximum scalar entropy S g (I\9), while A is the characteristic size parameter of the 
system. 
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Proof. Assuming that scalar entropy S g (I\9) and macroscopic observables I exhibits 
an extensive growth with the size parameter A, it is easy to check the following scaling 
dependencies: 

S 9 ~Al J ^~1/A,^~A, 

/~A J \ r*. ~ 1/A, R ljkl ~ 1/A 3 , R ~ 1/A. 1 ; 

Thus, the curvature scalar constitutes a size effect of order 1/A in the 

fluctuating behavior associated with the distribution function ([7]). According to 
Riemannian geometry [TJ, the asymptotic expression of the geometric rate a(r\6) for 
a small sphere can be expressed in terms of the curvature scalar R(I\9) at the center 
of the n-dimensional sphere (J, r) as follows: 

a(r\6) = 1 - — R(I\6)r 2 + 0(r 2 ), (86) 
6n 

which leads to the desirable result ( 1841) performing the integration of Eq.( l82l) . It 
is worth to remark that the estimation formula (|84p is applicable whenever the 
curvature scalar R{J\6) be sufficiently small. This observation suggests its possible 
applicability beyond the framework of the short-range interacting systems whenever 
the gaussian approximation (1321) turns a licit treatment to describe the statistical 
behavior of the macroscopic observables /. ■ 

5. Final remarks 

As already evidenced, Riemannian formulation of classical statistical mechanics 
arises as a suitable extension of Ruppeiner geometry of thermodynamics. A main 
consequence of this proposal is the interpretation of the scalar entropy S g (I\9) and 
its associated thermo-statistical quantities as purely geometric notions derived from 
the Riemannian structure of the manifold of macroscopic observables M.q. Besides, 
the non-vanishing of the statistical curvature tensor Rijki(I\6) of the manifold M. e 
constitutes a direct indicator about the existence of irreducible statistical dependence 
for the equilibrium distribution function dp(I\6) associated with the covariant 
redefinition of the Einstein postulate (J7J). 

At first glance, thermo-statistics geometry shares some analogies with Einstein 
gravitation theory. In particular, the notion of statistical correlations could be 
regarded as the statistical counterpart of the concept of interaction. Thus, the 
relation between statistical curvature tensor R%jki{I\0) and the existence of irreducible 
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statistical dependence for the equilibrium distribution function (J7J) is analogous to 
the relation between the spacetime curvature tensor Rijki(t,r) and the irreducible 
character of gravitation (while the existence of inertial forces can be avoided with 
an appropriate selection of the reference frame, gravitational forces are irreducible 
because of this universal interaction is just the manifestation of the spacetime 
curvature). 

Before to end this section, it is worth to comment some open questions that 
deserve a detailed analysis in forthcoming contributions. The present work inspires 
a geometric reformulation of classical fluctuation theory, which should lead to the 
derivation of a certain set of covariant fluctuation theorems and dynamical equations 
to describe the system relaxation. Similarly, it is worth to analyze the implications 
of thermo-statistics geometry on the study of phase transition phenomena as well 
as its possible extension to the framework of quantum distribution functions. 
According to the relationship between statistical curvature and the existence of 
irreducible statistical dependence for the distribution functions, an extension of the 
present proposal to quantum statistical distributions should provide an alternative 
characterization for the notion of quantum entanglement. 
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